A Measurement of sin 2 6*w in z/N Scattering from NuTeV 
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The NuTeV experiment at Fermilab presents a determination of the electroweak mixing angle. 
High purity, large statistics samples of v^N and V^N events allow the use of the Paschos-Wolfenstein 
relation. This considerably reduces systematic errors associated with charm production and other 
sources. With Standard Model assumptions, this measurement of sin 2 Ow indirectly determines the 
W boson mass to a precision comparable to direct measurements from high energy e + e~ and pp 
colliders. NuTeV measures sin 2 6> w (on ~ shcll) = 0.2253 ± 0.0019(stat.) ± 0.0010(syst.), which implies 
M w =80.26 ±0.11 GeV. 
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Neutrino scattering experiments have contributed to our understanding of electroweak physics for more than three 
Q^ ■ decades. Early determinations of sin 2 Ow served as the critical ingredient to the Standard Model's successful predic- 
tion of the W and Z boson masses. More precise investigations in the late 1980's set the first useful limits on the 
top quark mass. Just as early measurements contributed to the accurate predictions of the W, Z, and top quark 
masses before their direct observation, recent results fron 
pp collider data likewise constrain the Higgs boson mass. 



masses before their direct observation, recent results from neutrino-scattering experiments combined with e + e and 



The measurement presented here represents the most precise determination of the electroweak mixing angle from 
neutrino-nucleon scattering to date. The result is a factor of two more precise than the previous most accurate i/N 
measurement fjj. 

II. METHODOLOGY 

In deep inelastic neutrino-nucleon scattering, the weak mixing angle can be extracted from the ratio of neutral 
current (NC) to charged current (CC) total cross sections. Previous measurements relied on the Llewellyn- Smith 
formula, which relates these ratios to sin 2 Ow for neutrino scattering on isoscalar targets 0] : 
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The above equations are exact only for tree level scattering off an idealized isoscalar target composed of first generation 
light quarks. Corrections must be made for the non-isoscalar target, the heavy quark content of the nucleon, radiative 
effects, higher- twist processes, the longitudinal structure function {Rl), and charm production. This last effect is 
most important. Unfortunately, previous determinations of sin 8\y measured in this way were subject to the same 
charm production uncertainties (resulting from imprecise knowledge of the charm quark mass) that dominated the 
CCFR error [fjj. This ultimately limited the precision of neutrino measurements of electroweak parameters. 

An alternate method for determining sin 8w that is much less dependent on the details of charm production and 
other sources of model uncertainty is derived from the Paschos-Wolfenstein quantity, R~ M: 
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Because R~ is formed from the difference of neutrino and antincutrino cross sections, almost all sensitivity to the effects 
of sea quark scattering cancels. This reduces the error associated with heavy quark production by roughly a factor of 
eight relative to the previous analysis. The substantially reduced uncertainties, however, come at a price. Unlike R", 
the ratio R~ is more difficult to measure experimentally because neutral-current neutrino and antineutrino events 
have identical observed final states. The two samples can only be separated via a priori knowledge of the incoming 
neutrino beam type. 

III. THE NEUTRINO BEAM AND THE NUTEV DETECTOR 

High-purity neutrino and antineutrino beams were provided by the Sign Selected Quadrupole Train (SSQT) at the 
Fermilab Tevatron during the 1996-1997 fixed target run. Neutrinos are produced from the decay of pions and kaons 
resulting from interactions of 800 GeV protons in a BeO target. Dipole magnets immediately downstream of the 
proton target bend pions and kaons of specified charge in the direction of the NuTeV detector, while wrong-sign and 
neutral mesons are stopped in beam dumps. The resulting beam is almost purely neutrino or antineutrino depending 
on the selected sign of the parent mesons. The measured v^ contamination in the v^ beam is less than 1/1000 and 
the v^ contamination in the V ^ beam is less than 1/500. In addition, the beam is almost purely muon- neutrino with 
a small contamination of electron neutrinos (1.3% in neutrino mode and 1.1% in antineutrino mode). 

Neutrino interactions are then observed in the NuTeV detector, which is located approximately 1.5 km downstream 
of the proton target. The detector consists of an 18m long, 690 ton steel-scintillator target followed by an instrumented 
iron-toroid spectrometer (Figure p]). The target calorimeter is composed of 168 3 m x 3 m x 5.1 cm steel plates 
interspersed with liquid scintillation counters (spaced every two plates) and drift chambers (spaced every four plates) . 
The scintillation counters provide triggering information as well as a determination of the longitudinal event vertex, 
event length and visible energy deposition. The mean position of hits in the drift chambers help establish the transverse 
event vertex. The toroidal spectrometer, which determines muon sign and momentum, is not directly used for this 
analysis. 

LAB-E Detector - Fermilab E815 (NuTeV) 
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FIG. 1. The NuTev detector. 



The detector was continuously calibrated through exposure to a wide energy range of test beam hadrons, electrons 
and muons delivered in a separate beamline during each accelerator cycle. 



IV. EVENT SELECTION 

For inclusion in this analysis, events must deposit at least 20 GeV of visible energy in the detector. This cut ensures 
full efficiency of the event trigger, proper vertex determination, and reduction of cosmic ray background events. Events 
must also have a vertex that lies within a 1.0 m box around the center of the detector which is at least 0.4 m of 
steel from the upstream end and 2.4 m of steel from the downstream end of the detector. This "box cut" ensures 
shower containment and minimizes the contribution from v e interactions. The chosen longitudinal fiducial volume 
ensures both that the event is neutrino-induced and that a meaningful event length can be measured. After all cuts, 
the total analysis sample consists of 1.3 million neutrino and 0.3 million antineutrino events with a mean energy of 
approximately 125 GeV. 

V. EVENT SEPARATION 

In order to measure sin 8w, observed neutrino events must be separated into charged current (CC) and neutral 
current (NC) categories. Figure shows typical CC and NC candidate events in the NuTeV detector. 




FIG. 2. Event displays of a NuTeV charged current candidate event (top) and neutral current candidate (bottom). 



Both CC and NC neutrino interactions initiate a cascade of hadrons in the target that is registered in both the 
scintillation counters and drift chambers. Muon-neutrino CC events are distinguished by the presence of a final 
state muon. The muon typically penetrates well beyond the hadronic shower and deposits energy characteristic of a 
minimum- ionizing particle in a large number of scintillation counters. The muon track is clearly visible in the CC 
event at the top of Figure 0. On the other hand, there is no final state muon in muon-neutrino NC events; the final 
state neutrino is invisible, so v^ NC events can only be recognized by their hadronic shower. 

Given the differing event topologies, one can statistically separate CC and NC interactions on the basis of event 
length (i.e., on the presence or absence of a muon in an event). The length of an event is based on the longitudinal 
energy deposition in the calorimeter and is simply the number of scintillation counters spanned by the event. Events 
with a long length (spanning > 20 counters) are identified as CC candidates and events with a short length (spanning 
< 20 counters) are identified as NC candidates. The experimental quantity that is measured in both modes is the 
ratio: 

# SHORT events # NC candidates 
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# LONG events # CC candidates 

Neither sample is of course entirely pure. Long events are predominantly v^ CC interactions but contain small 
contaminations from both v^ NC events where the hadronic shower fluctuates to long lengths, as well as muons which 



are produced in upstream neutrino interactions in the shielding and later undergo a hard bremsstrahlung interaction 
in the detector. Short events are primarily v^ NC interactions but also contain contributions from low energy v^ CC 
events in which the muon ranges out, wide angle v^ CC events in which the muon exits the side of the detector, v e 
CC events and cosmic rays. The ratios (R mea s) of short to long events measured from the NuTeV data, are 0.4198 ± 
0.0008 in the neutrino beam and 0.4215 ± 0.0017 in the antineutrino beam. A Standard Model value of sin 2 By/ can 
be directly extracted from these measured ratios by using a detailed Monte Carlo simulation of the experiment. The 
Monte Carlo must include the integrated neutrino fluxes, the neutrino cross section as well as a detailed description 
of the NuTeV detector. The following sections will discuss these three components. 

VI. THE MUON-NEUTRINO AND ELECTRON-NEUTRINO FLUXES 

The Monte Carlo simulation requires neutrino flux information as input. In particular, a precise determination of 
the electron neutrino contamination in the beam is essential. The measured short to long ratios are directly impacted 
by the presence of electron neutrinos in the data sample because v & charged current interactions, which usually lack 
an energetic muon in the final state, are almost always identified as neutral current interactions in the detector. 

The electron neutrino flux is estimated using a detailed beam Monte Carlo. Most of the observed v e 's result from 
K^ — * 7r°e ± i/ e decay. The beam simulation can be tuned to describe v e and v e production from charged kaon decay 
with high accuracy given that the K decay contribution is well constrained from measurements of the observed 
Vp and v ^ fluxes (as shown in Figure pi). Because of the precise alignment of the beamlinc elements and the low 
acceptance for neutral particle propagation in the SSQT, the largest uncertainty in the calculated electron neutrino 
flux results from the 1.5% uncertainty in the K^ —* 7r°e ± ^ e branching ratio. The result is a factor of three reduction 
in the electron-neutrino flux uncertainty when compared to CCFR. 
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FIG. 3. v^ and v^ energy spectra comparing data (points) and the tuned beam Monte Carlo (histogram) 



VII. CROSS SECTION MODEL 



Neutrino-nucleon deep-inelastic scattering processes are simulated using a leading order cross section model. The 
cross section model incorporates leading order parton momentum distributions measured using the same target and 
cross section model as NuTeV Q. Small modifications adjust the parton densities to produce the inherent up-down 
quark asymmetry consistent with muon scattering |5|] and Drell-Yan J6[ data. A leading order analysis of dimuon 
events in CCFR |7J] provides the shape and magnitude of the strange sea. Mass suppression from charged-current 
charm production is modeled using a leading order slow-rescaling formalism whose parameters are measured from the 
same high-statistics dimuon sample. Electroweak and QED radiative corrections to the scattering cross sections are 
applied using radiative correction routines supplied by D. Yu. Bardin 0. Higher twist contributions are quantified 
using a Vector-Dominance-Model as implemented by Pumplin |13[ and constrained by lepto-production data [[12| . A 
global analysis performed by L. Whitlow pj] provides a parameterization of the longitudinal structure function, R^. 



VIII. DETECTOR RESPONSE 

The Monte Carlo must also accurately simulate the response of the detector to the product of neutrino interactions in 
the target. The critical detector parameters that must be modeled are the calorimeter response to muons, measurement 
of the neutrino interaction vertex, and the range of hadronic showers in the calorimeter. Precise determination of the 
various detector effects are made possible through extensive use of both neutrino data and large samples of calibration 
beam data. The efficiency, noise, and effective size of the scintillation counters are measured using neutrino data or 
test beam muons. Vertex finding resolutions and biases are studied using neutrino data combined with information 
from a detailed GEANT-based simulation of the detector. Test beam pion data provides information on hadronic 
shower lengths. 

IX. EVENT LENGTH DISTRIBUTIONS 

An important test of the Monte Carlo is its ability to predict the length distribution of events in the detector. 
Figure shows event length distributions of the final data sample compared to the Monte Carlo prediction assuming 
the best- fit sin 2 W (see next section). Events reaching the toroid, which comprise about 80% of the CC sample, have 
been left out for clarity but are included in the normalization of the Monte Carlo to the data. Data and MC agree 
well in each beam mode. 
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FIG. 4. Event length comparisons for neutrino and antineutrino events. The solid curve is the Monte Carlo prediction. The 
neutral current to charged current event separation (indicated by the arrow) is made at a length of 20 counters, approximately 
2m of steel. 



X. EXTRACTION OF sin 2 6> w 

Given separate high-purity neutrino and antineutrino data sets, NuTeV measures the following linear combination 
of R v and R v : 



R~ =R V -xR", 



(0) 



where x is chosen to be 0.5316. This value is obtained using Monte Carlo to minimize the dependence of R~ on 
the charm quark mass. This approach explicitly minimizes uncertainties related to the suppression of charm quark 
production, largely eliminates uncertainties related to sea quark scattering, and reduces many of the theoretical and 
detector uncertainties common to both the neutrino and antineutrino samples. The single remaining free parameter 
in the Monte Carlo, sin 2 9yy, is then varied until the model calculation of R~ agrees with what is measured in the 
data. The preliminary result from the NuTeV data sample for Mt op =175 GeV and Afniggs=150 GeV is: 
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0.2253 ± 0.0019(stat.) ± 0.0010(syst.) 
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Leading terms in the one-loop electroweak radiative corrections ||] to the W and Z self-energies produce a small 
residual dependence of our result on M top and Mjuggs- We explicitly quote the dependence of our measurement on 
the top and Higgs masses. The effect is small given the existing uncertainty on the top quark mass and the logarithmic 
dependence on the Higgs mass. For example, the total change in sin 2 9\y from varying the Higgs mass from 100 GeV 
to 1 TeV is only 0.0005 (as is indicated by the shaded gray bands outlining the NuTeV result in Figure Wj. 
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Having chosen the convention, sin 2 #^ n = 1 — -g%-, where Mw and Mz are the physical gauge boson masses, 

our result implies: 

M w = 80.26 ± 0.10(stat.) ± 0.05(syst.) = 80.26 ± O.llGeV (9) 

A comparison of this result with direct measurements of Mw is shown in Figure |5j. Our measurement is in good 
agreement with Standard Model expectations and is consistent with the most recent measurements from W and Z 
production as well as from other neutrino experiments (Figure |6|) . 
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FIG. 6. Experimental constraints presented on 
the Mw-Mtop plane (at the time of DPF99). The 
two narrow lines indicate the Standard Model pre- 
dictions for M H i S ss = 100 and 1000 GeV. 



If these plots are now updated to include new measurements from CDF ||, D0 [JHJ and LEP |ll|] that were 
presented at Electroweak Moriond in March of this year, it is interesting to note that tighter constraints on the Higgs 
mass are beginning to emerge (Figure @) . 
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FIG. 7. Direct W boson mass measurements 
compared with this result (updated to include new 
measurements presented at EW Moriond99). 
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FIG. 8. Experimental constraints presented on 
the Mw-Mtop plane (updated to include new mea- 
surements presented at EW Moriond99). 



XI. CONCLUSIONS 

NuTeV has successfully completed its data-taking and has extracted a preliminary value of sin Ow- The precision of 
this result is a factor of two improvement over previous measurements in vN scattering because of reduced systematics 
associated with measuring the Paschos Wolfenstein ratio, R~ . Interpreted within the framework of the Standard 
Model, this result is equivalent to a determination of the W mass and is consistent with direct measurements of Mw- 
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